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Q Abstract. For a certain class of partitions, a simple qualitative relation is 

^~^ observed between the shape of the Young diagram and the pattern of zeroes of 

^^ the Wronskian of the corresponding Hermite polynomials. In the case of two- 

^^ term Wronskian W{Hn, ii^n+fc) we give an explicit formula for the asymptotic 

^S) shape of the zero set as n — > oo . Some empirical asymptotic formulas are given 

for the zero sets of three and four-term Wronskians. 
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1. Introduction 
Consider a Schrodinger operator 



d 



•2 



with a rational potential u(^z), not necessarily real. Such an operator is called 

^ monodromy-free if all the solutions of the corresponding Schrodinger equation Lip — 

lO Xtp are meromorphic in the whole complex plane for all A. 

^^ The first classification result here is due to Duistermaat and Griinbaum |2], 

VO .... 

jT^ who described all monodromy-free operators with rational potentials decaying at 

infinity. The pole configurations of the corresponding potentials were studied earlier 
)^ by Airault, McKean and Moser 1 in relation to rational solutions of the KdV 

equation. 

Oblomkov [31 generalised Duistermaat-Griinbaum's result to the quadratic growth 
case. He showed that all the rational monodromy-free operators with rational po- 
tentials growing as z^ at infinity are the results of Darboux transformations applied 
to the harmonic oscillator. The corresponding potentials have the form 



S .,,.^_ .d' 



uiz) = -2— log WiHk, ,...,HkJ + z'+c, 

where Hk{z) is the fc-th Hermite polynomial, /ci > ^2 > • • • > fc„ is a sequence of dif- 
ferent positive integers and W(/i, . . . , /„) is the Wronskian of functions /i, . . . , /„. 
We are interested in the geometry of the pole configurations of the corresponding 
potentials {locus in the terminology of Airault, McKean and Moser), which are the 
same as the zero sets of the corresponding Wronskians. This locus has an interesting 
relationship with the Calogero-Moser problem and log-gas in a harmonic field, see 
[3] . In the case when /ci , . . . , fc„ are consequent numbers it can be also interpreted 
as the pole set of some rational solutions of the fourth Painleve equation and has a 
regular rectangle-like structure in the complex plane, as was revealed numerically 
by Clarkson [S] . A natural question is what kind of pattern do we have in general. 
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Let US label these potentials by partitions A = (Ai, . . . , A„), Ai > • • • > A„ > 1, 
such that Xi = ki — n + i, i — 1, . . . ,n : 

ki = Xi + n - 1, k2 = X2 + n - 2, . . . , fc„_i = A„_i + 1, fc„ = A„. 

Our main observation (based on numerical experiments using Mathematica) is 
that although for a general partition A the picture is quite complicated, for its 
doubled version 

A2x2 = ((2Ai)2,...,(2A„)2) 

there exists a simple qualitative relation between the shape of the Young diagram 
and the pattern of zeroes of the corresponding Wronskian. 

In the case of the two-term Wronskian W{Hn, Hn+k) we have some quantitative 
results. Namely, for fixed k and large n we give an explicit formula for the curve 
the scaled zeroes w = z/^/2n = u + iv in the region |u| < 1 — ^, \v\ > £-^^^^, e,5 > 
lie on asymptotically : 

(1) H = ^^^^(ln(f)+ln(l-.^) + lln|l-T|(.)|), 

where T/j (x) = cos k arccos x is the fc-th Chebyshev polynomial. The distribution 
of the real parts of the zeroes on this curve satisfy Wigner's semicircle law. The 
derivation is based on a version of the classical Plancherel-Rotach formula J6 found 
by Deift et al in [7]. We give also some empirical formulas for the three and four- 
term Wronskians. 

2. Wronskians of Hermite polynomials and their zeroes 

Hermite polynomials Hn{x) are the classical orthogonal polynomials with Gauss- 
ian weight w{x) = e~^ (see e.g. [B]). They can be given by the formula 

and satisfy the recurrence relation 

Hn+i{x) = 2xHn{x) - 2nHn^i{x). 
Here are the first few of them: 

Ha{x) = I, Hi{x) = 2x, H2{x) = Ax^ - 2, Hi{x) = 8a;^ - 12a;, 

H^{x) = 16a;'' - A^x^ + 12, H^{x) = 32a;^ - leOa;^ + 120a;, ... 
We are using the normalisation where the highest coefficient of i/„ is 2", but this will 
not be essential in what follows. What is important for us is that ipn = Hn{x)e^^ /^ 
are the eigenfunctions of the harmonic oscillator: 

(-^ + x^)^Im = {n + 1/2)Vj„, , n = 0, 1, . . . . 
Let A = (Al, . . . , A„) be a partition and consider the Wronskian 

The Wronskians W\ have the following properties: 

1. W\{z) is a polynomial in z of degree |A| = Ai + A2 + • • • + A„, 

2. Wx{-z) = (-I)I^IW^a(^), 

3. Wx>{z) — (— i)I^WA(iz), where X* is the conjugate of X. 
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Figure 1. Ferrers diagram for the partitfon A = (5,3,3,1). Left: 
standard version. Right: French version 



Figure 2. Standard and French versions of the Young diagram 
for the partition A — (5, 3, 3, 1). 



Recall that the conjugate to a partition A is a new partition, whose Young di- 
agram is the transpose of the diagram of A. To prove the last (duality) property 
we note that the harmonic oscillator has also the following (growing, hence formal) 
eigenfunctions V^ = H^{x)e^ /^, H^{x) — (— i)"iJ„(ix) with the negative eigen- 
values A = —n + 1/2, n = 0, 1, . . . . The claim is that the result of applications of 
Darboux transformations at the levels Xi +p — i, i = 1, . . . ,p, p is the length of the 
partition, and at the negative levels j — q — X*, j = 1, ... ,q,q is the length of the 
conjugate partition, lead to the potentials different only by a constant shift. This 
follows, for example, from proposition (1.7) in Macdonald's book (8j), or from the 
so-called Maya representation of the partition. 

Recall the following well-known diagrammatic representations of a partition A 
(see e.g.[H]). The first one is the set of points (i, j) G Z^ such that 1 < j < A^. 
Following [5] we will call it a Ferrers diagram. There are two ways to draw this. 
One convention (motivated by matrix theory) is that the row index i increases 
downwards while j increases as one goes from left to right. Another way (sometimes 
called French) is to use a natural Cartesian coordinate representation (see Fig. IT]). 

The second, most common, way, known as a Young diagram, is to use boxes 
rather than bullets (see Fig. [2|. 

Since the Wronskians are labelled by the partitions, we can ask a natural ques- 
tion: how is the geometry of the corresponding diagram of A related to the pattern 
of the zeroes of Wx{z) ? Figure [3J produced with the help of Mathematica, shows 
that in general such a relation is not easy to see. Another example is the partition 
A = (n, n — 1, n — 2, . . . , 2, 1) with a triangular Young diagram, for which the cor- 
responding Wronskian Wx (up to a multiple) is simply 2:"("+i)/2, so we just have 
one zero at z = with multiplicity n{n + l)/2. 
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Figure 3. Zeroes of the Wronskian Wx with A = (28, 16, 10, 6, 4, 4, 3, 1). 



Figure 4. Diagram of the doubled partition A^^^ for A = (5, 3, 2). 



This is why we found it very interesting that for a special class of partitions, 
which we call doubled, one can read off the partition from the pattern of zeroes in 
a straightforward way. 



3. Doubled partitions and their diagrams 

Let A = (Ai, A2, . . . , A„) be a partition. Define its doubled version as 

A = (2Ai, 2Ai, 2A2, 2A2, . . . , 2A„, 2A„). 

In other words, we double all parts and take them twice. For example, when 
A = (5,3,2) the doubled version is A^^^^ ^ (10,10,6,6,4,4) = (10^,62,4^), where 
the power denotes how many times this part is repeated. 

Note that the shape of the Young diagram of the doubled version is similar 
(with scaling factor 2) to the initial one. However, for the doubled partitions there 
is another natural way to represent them, which combines both the usual and the 
French ways. Namely one can put the diagram of A in all 4 quadrants as in Fig. 
[4] One can naturally define the Ferrers version, which we combine with the Young 
version by putting bullets at the centre of each box. 
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Figure 5. BuUeted diagram of the doubled partition A^^^ for A 
(5, 3, 2) and the zeroes of the corresponding Wronskian W\2x2 



J.O- 



FiGURE 6. The same comparison for A = (7,4, 1). 



2x2 



Our main observation is that the diagram of the doubled partition X^^^ gives a 
good qualitatative description of the zero set of the corresponding Wronskian W\2x2, 
see Fig. [5]and|6] We beheve that this works for any partition A with distinct A^. 

When some parts are equal then we may have interference between the rows of 
corresponding zeroes, see Fig. [7j 



• 1 - • 

• • • • 

• • • • 

-2-1 1 2 

• • • • 



Figure 7. Interference between the rows of zeroes of W\2x 
partition A = (5, 3, 3, 1) with two equal parts. 



for 



One can generalise this relation to the case of half-integer partitions (Ai, . . . , A„) 
with some of the parts being half-integers. An example is A = (1,5/2,11/2) for 
which the doubled partition is A^^^ = (2^, 5^, 11^). The corresponding Ferrers dia- 
gram has some bullets on the vertical axis and gives a good qualitative picture of 
the zero set of the corresponding Wronskian (see Fig. [s]). 
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Figure 8. Ferrers diagram (left) of A^^^ and zeroes of W\2x2 
(right) for half-integer partition A = (11/2,5/2, 1). 

The following asymptotic analysis shows the limits of this comparison already 
in the simplest case of one-row Young diagram A. 



4. Asymptotic behaviour of zeroes of two-term Wronskians 

Consider now the two-term Wronskian W{Hn, Hn+k), corresponding to the par- 
tition A = {n + k—l,n),k> 1. Let us fix k and let n — > oo. To study this behaviour 
of zeroes in this limit we can use the following version of Plancherel-Rotach formula 
due to Deift et al [7]. [] 

In the scaled variable w = z/\/2n there are several regions with different asymp- 
totic behaviour of the Hermite polynomials (see Fig. 9). The most relevant for us 
is the region Bs, where we have the following asymptotics 



A, 




B, 




-l-S 



-1 + 6 



1-5 



Figure 9. Asymptotic regions in scaled variable w 



H, 



2 / 1 1 

(z)e"^ ^ C„(l-u;2)-i ( cos{2n<d{w)+x{w)){^+0{-))+sm{2n<d{w)~x{w))0{-) 
\ n n 



with Cn — \/ ^{2n) ■>■ , 9(w) — ^wVl — w'^ + ^ arcsinw— ^ and x(w) — \ arcsinu; 

(see [71). 

Using this we can show that 

W^(ff„,i/„+fe) = -^ 



sin A I. 



(sin Afe + k sin(2$ + Afe))! (1 + 0(-)), 

J n 



4n(l - w2) 

where $ = 2n0 -|- X: ^kiw) = 2kQ(w) — kwQ'{w) = fcarccosw. In the upper 
half-plane we have two competing terms: sin A^. term and the negative exponential 



We are very grateful to Ken McLaughlin for attracting our attention to this important paper 
during " Dubrovin-60" conference at Sardinia in June 2010. 
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component e^*(^*+^'''' of sin(2<i> + A^). For w = u + iv with small v « 1 we can 
approximate $(u)) as 

$(u + iv) w 2nQ{u) + 2nivQ' [u) = 2nQ{u) + 2niv\/l - v? 

since 0'(w) = v 1 — w"^ . Equating the moduli of the two competing terms, we have 

k 



sinAi-fwl 



8n(l 



_ Anv^/l — u^ 



(2) V = (ln(-^) +ln(l -u^) + In I sin(fc arccosu)!). 

4nVl -u^^ k ' J 

This leads to the formula (IT]) for the curve on which the scaled zeroes lie asymp- 
totically as n — >■ cx) in the region |u| < 1 — (5, \v\ > e-^^^Jl^ ^^§ > q. Comparing 
the arguments of the leading terms and using the relation dQ = y/l — w'^dw we 
see that the real parts of the corresponding zeroes are distributed according to the 
famous Wigner's semicircle law from the random matrix theory |12j : the number 
Na,l3{n, k) of scaled zeroes w — u-\-iv oi W{Hm Hn+k) hi the upper- half plane with 
the real parts in the interval (a,/3) satisfies 

(3) hm ^^-lil^ . ^ r VY^Uu. 

n-i-oo n TT J^ 

When w = we have fc — 1 real (scaled) zeroes Um = x„il \/2n asymptotically 
given by 1 — T^iu) ~ 0: as rt — > cx) 

(4) w,„ = — z= ^- cos — — , m = l,...,fc- 1. 

\/2n k 



The unsealed zeroes z — x + iy in the region 

n,,S : \x\ < (l ~ S)V2n, \y\ > 



logn 

lie asymptotically on the curve 

(5) M = — ;^l=(ln(^)+ln(l-f^) + iln|l-T|(^; 

2\/2n - x^^ k ' 2n 2 ^2n 

where Tk{x) is the k-th Chebyshev polynomial. Figure 10 shows a good agreement 
of this formula (curve) and numerical Mathematica calculation of zeroes (dots) in 
the case when n = 100, fc = 5. The 4 real zeroes in this case approximately are 

x « ±V200 COS -— = 5V2 -^^, m = l,2 

5 2 

in agreement with the picture. 

The case n = 100, k = I corresponds to the doubled partition (100, 100) = 50^^^. 
Figure 11 shows that shapes of the Young diagram and the corresponding zero set 
coincide only qualitatively. Indeed, the corresponding asymptotic curve in this case 
is not just two straight lines but given by ([5| with fc = 1: 
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Figure 10. Comparison in the case n=100, k=5 
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Figure 11. Zeroes of W^(iJioo,i?ioi) 



5. Empirical asymptotic formulas for 3-term and 4-term Wronskians 

5.1. Three-term Wronskians. One can prove the following identity for 3-term 
Wronskians of the eigenfunctions of harmonic oscillator 

We will restrict ourselves with the case k = I. In that case we have 

(6) W3 = k{tpn-kW{lpn, i^n+k) - Ipn+kW (tpn-k, -^n)) 

Using this, previous results about 2-term Wronskians and experimenting with 
Mathematica we can suggest the following empirical formula for the limiting shape 
of the non-real zeroes for large n and k « n: 

. . 1 



\V\ 



.(\n{^)+\n{l-x^/2n) + \\n\l-Tl{x/^/^)\) 



(7) 



■\/2n — x^ 

Mathematica plots of zeroes against the corresponding curve below show a good 
agreement with this formula. 
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Figure 12. Comparison of the zeroes and the curve for 
Left:P^(Hino, i?ioi, ^102) Right:W^(Hioo, -ffios, ^loe)- 




Figure 13. Comparison of the zeroes and the curve for 
W^(^ioo,-H^i05,^iio) ■ 



5.2. 4-term Wronskians. For Wronskians WA = M^('0„, V'n+fc, i^n+k+ui^n+k+i+m) 
one can show that 



(8) 



H/4 



l{k + l + m)W{iln, i>n+k)W {iln+k+hi'n+k+l+Tn) 
-kmW{lpn,i^n+k+l+ni)W{'lpn+k,'^Pn+k+l) 

Assuming first that k — I — m we have 
(9) WA = 3fc2w^(V„,V^„+fe)VF(V^„+2fc,^„+3fc) 

Mathematica plots suggest that the zeroes for large n and small k « n asymp- 
totically lie on two curves, for which we have the following empirical formulas: 



(10) 



H = -^=L==(ln(^)+ln(l-xV2n) + iln|l-T2(x/y2^)|) 



2^2^- a; 
for the middle curve and 
3 



(11) H^ ^^___ (ln(^)+ln(l-:rV2n) + ^ln|l-T,^(x/^/2;^)|) 

for the outside curve. 
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Below we compare these curves and Mathematica plots of zeroes for n — 100 and 
values of k ranging from 1 to 4. The first picture corresponds to the double partition 
A^^^ with A ~ (50, 50). In all cases we expect the real parts of properly scaled zeroes 
to satisfy Wigner's semicircle law. We note a peculiar behaviour of zeroes near the 
points where 1 — T^{x/\/2n) = 0, which needs a special investigation. 






Figure 14. Comparison of the zeroes and the curve for 
Leit:W{HiQQ, iJioi, i?io2, -ffios) Right :W^(iyioo, ^102, -^104, ^loe)- 








Figure 15. Comparison of the zeroes and the curve for 
Leit:W{Hioo, Hio3, HioQ, Hiog) Right:W(_ffioo,ffio4, ^1087^112)- 

Consider now the case k = m = 1 and W4: = W^(iJ„,_ff„+i,7?„+;+i, iJ„+/+2), 
corresponding to the doubled partitions A^^^ with A = (n/2, {I — l)/2). The empir- 
ical formulas for the asymptotic zero curves for large n are 



(12) 



1 



m 



2y/2n - x' 



:(ln4n+ -ln{l-x^/2n)- - \n \1 - T^ {x / V2n)\) 
V 2 2 / 



for the middle curve and 
, , 1 



(13) 



12/1 = 



V2r 



In 



^ In (1 - x'/2n) + J In |1 - T^{x/V2^)\) 



m - x'^ '- 51 2 I 

for the outside curve. They seem to work fairly well when I < n/A but not too small. 
The cases with n — 50, / = 11 and n = 60, / = 10 are shown below. We should say 
that at the moment this part is just experimental mathematics and requires further 
investigation. 
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Figure 16. Comparison of the zeroes and the curves for W{Hqq, H^i, Hq2, Hq^). 




Figure 17. Comparison of the zeroes and the curves for W{Hqo, Hqi, ifyi, -^72). 



6. Some conjectures 

The following property of the Wronskians of Hermite polynomials was conjec- 
tured by the third author in 1990s in relation with the corresponding locus problem 
solved by Oblomkov [3] . If this property holds it would give a way to a more effective 
proof of his result, which is still very desirable. 

Conjecture 1. For every partition X, all the zeroes of W\{z) are simple except 
possibly for z — 0. 

Note that the multiplicity m of 2; = for W\ can be easily computed and has 
the form 

d(d+l) 
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where d = p — q is the difference between the numbers p and q of odd and even 
elements respectively among the sequence Ai +n— f , A2 + n — 2, . . . , A„_i + 1, A„. In 
particular, for the triangular Young diagram A = (n, n — 1, ... ,2,1) we have d = n 
and m = n(n + l)/2 = degWx, so the corresponding Wronskian W\ = C„2:"("+^'/^ 
and all the zeroes collide at zero. 

An interesting question is if the number of real zeroes of Wx{z) can be effectively 
described in terms of the corresponding Young diagram. For doubled partitions we 
have the following conjecture. 

Conjecture 2. For doubled partitions v — (/i^, . . . , /i^) with distinct parts, the 
Wronskian W^{z) has no real roots and has as many pure imaginary roots as there 
are odd numbers among fii, . . . , fin- 

In the special case when n — 1 and v = (m, m) we can prove this using the 
integral representation of the corresponding Wronskian known from the random 
matrix theory (see Brezin-Hikami [TO]): 

-°° "'-°° i<j fc=i 

Finally, it would be very interesting to understand how special the Hermite 
polynomials are and how much of this can be generalised to other orthogonal poly- 
nomials and to the sextic growth case [IT] . 

7. Acknowledgments 

One of us (APV) is grateful to the Institute for Mathematical Research at ETH 
Zurich for the hospitality in April 2010 and to Robert Milson for stimulating dis- 
cussions. 



References 

H. Airault, H.P. McKean, J. Moser Rational and elliptic solutions of the Korteweg-de Vries 
equation and a related many-body problem. Comm. Pure Appl. Math. 30 (1977), 95-148. 
J.J. Duistermaat and F.A. Griinbaum Differential equations in the spectral parameter. Comm. 
Math. Phys., 103 (1986), 177-240. 

A. A. Oblomkov Monodromy-free Schrodinger operators with quadratically increasing poten- 
tials. Theor. Math. Phys., 121 (1999), 1574-1584. 

A. P. Vcsclov On Stieltjes relations, Painleve-IV hierarchy and complex monodromy. J. Phys. 
A 34 (2001), no. 16, 3511-3519. 

P.A. Clarkson The fourth Painleve equation and associated special polynomials. J. Math. 
Phys. 44 (2003), 5350-5374. 

G. Szego Orthogonal polynomials, American Mathematical Society Colloquium Publications, 
Vol XXIII, 1959. 

P. Deift, T. Kriecherbauer, K. McLaughlin, S. Venakides, X. Zhou Strong asymptotics of 
orthogonal polynomials with respect to exponential weights. Comm. Pure Appl. Math., Vol. 
LII (1999), 1491-1552. 

I. Macdonald Symmetric functions and Hall polynomials. 2nd edition, Oxford Univ. Press, 
1995. 

B.E. Sagan The symmetric group. Representations, combinatorial algorithms, and symmetric 
functions. 2nd edition. Graduate Texts in Mathematics, 203. Springer- Verlag, New York, 
2001. 

E. Brezin, S. Hikami Characteristic polynomials of random matrices. Comm. Math. Phys. 
214 (2000), no. 1, 111-135. 

J. Gibbons, A. P. Veselov On the rational monodromy-free potentials with sextic growth. J. 
Math. Phys. 50 (2009), no. 1, 013513, 25 pp. 



ZEROES OF WRONSKIANS OF HERMITE POLYNOMIALS AND YOUNG DIAGRAMS 13 

[12] E. Wigner On the distribution of the roots of certain symmetric matrices, Ann. of Math. 67 
(1958), 325-328. 

Department of Mathematics, ETH Zurich, 8092 Zurich, Switzerland 
E-mail address: f elderSmath.ethz.ch 

Department of Mathematical Sciences, Loughborough University, Loughborough, 
Leicestershire, LEll 3TU, UK 

E-mail address: A. D. Hemery01boro.ac.uk 

Department of Mathematical Sciences, Loughborough University, Loughborough, 
Leicestershire, LEll 3TU, UK and Moscow State University, Moscow 119899, Russia 
E-mail address: A.P.VeselovOlboro.ac.uk 



